Boundary effects on the local density of states of one-dimensional 
Mott insulators and charge density wave states 



Dirk Schuricht, 1, 2 Fabian H. L. Essler, 3 Akbar Jaefari, 4 and Eduardo Fradkin 4 



'institute for Theory of Statistical Physics, RWTH Aachen, 52056 Aachen, Germany 
2 JAR A- Fundamentals of Future Information Technology 
The Rudolf Peierls Centre for Theoretical Physics, 
University of Oxford, 1 Keble Road, 0X1 SNP, Oxford, UK 
^Department of Physics, University of Illinois at Urbana- Champaign, 
1110 W. Green St, Urbana, Illinois 61801-3080, USA 
(Dated: November 16, 2011) 



We determine the local density of states (LDOS) for spin-gapped one-dimensional charge density 
wave (CDW) states and Mott insulators in the presence of a hard-wall boundary. We calculate 
the boundary contribution to the single-particle Green function in the low-energy limit using field 
theory techniques and analyze it in terms of its Fourier transform in both time and space. The 
boundary LDOS in the CDW case exhibits a singularity at momentum 2/cf, which is indicative of 
the pinning of the CDW order at the impurity. We further observe several dispersing features at 
frequencies above the spin gap, which provide a characteristic signature of spin-charge separation. 
This demonstrates that the boundary LDOS can be used to infer properties of the underlying bulk 
system. In presence of a boundary magnetic field mid-gap states localized at the boundary emerge. 
We investigate the signature of such bound states in the LDOS. We discuss implications of our 
results on STM experiments on quasi-lD systems such as two-leg ladder materials like Sri4Cu2404i. 
By exchanging the roles of charge and spin sectors, all our results directly carry over to the case of 
one-dimensional Mott insulators. 



Scanning tunneling microscopy (STM) and spectroscopy (STS) methods have proved to be a useful tool for studying 
strongly correlated electron systems such as carbon nanotubesi, high temperature superconductors (HTSC)2~— and 
rare-earth compounds^. STM experiments measure the tunneling current / between the sample and the STM tip as 
a function of its position x and the applied voltage V. This current can be expressed in terms of the local densities 
of states (LDOS) in the sample N(E,x) and the tip N tip (E) as£ 



where f(E) denotes the Fermi function. Assuming a structureless density of states in the tip, iVtip = const, this gives 
the following expression for the local tunneling conductance 



Eq. ([5]) shows that the tunneling conductance is proportional to the thermally smeared N(E, x) of the sample at 
the position of the tip. A non-trivial spatial dependence of the LDOS arises in presence of impurities. These break 
translational invariance and lead to a modification of the LDOS in their vicinity, from which one can infer characteristic 
properties of the bulk state of matter as well as the nature of its electronic excitations. Spatial modulations of the 
LDOS can be analyzed in terms the Fourier transform of the tunneling conductance. It follows from (J2J that this 
quantity is directly proportional to the corresponding Fourier transform of the LDOS, N(E,Q). This method of 
analyzing STM data was used very successfully^ to study quasiparticle interference in Bi 2 Sr 2 CaCu208+5. The spin 
dependence of the LDOS has also been investigated using magnetic tipsi Theoretical studies of STS have focused 
in particular on Luttingcr liquids^ and HTSC— r— . In the Luttinger liquid case an impurity has the same effects at 
low energies as a physical boundary^, which motivated studies of the LDOS in the vicinity of a chain end. The case 
of strongly correlated one-dimensional (ID) systems with spin or charge gaps is of considerable interest as well and 
pertains to quasi ID charge density wave (CDW) systems^ and Mott insulator a 14 ' 15 , carbon nanotube*^, (doped) 
two-leg ladder material a 17 ' 18 , and the stripe phases of HTSCiS. As compared to the Luttinger liquid case the presence 
of an interaction-induced gap makes these problems much more difficult to treat theoretically. In the following we will 
determine the LDOS for the low-energy limit of ID CDW states and Mott insulators in presence of a single boundary. 
The latter can be thought of as arising as the result of the presence of a strong potential impurity. Alternatively 
one can imagine inducing a boundary in a two tip STM setup, where the first tip is used to induce a boundary by 



I. INTRODUCTION 
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applying a high voltage and the LDOS is then measured with a second tip. A short summary of our results has 
appeared previously^. 

The outline of this paper is as follows: in section Ql] we present the field theory limit of ID CDW states and Mott 
insulators in presence of a single boundary. In section IIIII we summarize our results for the single-particle Green 
function. These are then used to determine the Fourier transform N a (E, Q) of the local density of states for hard- wall 
boundary conditions in section IPV1 Signatures of spin and charge excitations visible in N a (E,Q) are discussed in 
some detail. The effects of more general boundary conditions, including the formation of boundary bound states, 
are described in section [V] Section I VII deals with the effects of finite temperatures and implications of our results 
for STM experiments are discussed in section IVII1 The technical details of our calculations are presented in several 
appendices. 



II. THE MODEL 



Our analysis of the LDOS is based on the continuum description of certain ID CDW states and Mott insulators. 
The resulting quantum field theory for both cases is known as the U(l) Thirring models (with two "flavors".) The 
latter is known to arise as the effective low-energy description of a number of lattice models of spin-1/2 electrons as 
we discuss next. 

1. Half- filled repulsive Hubbard model 2 ^: 

This is the standard model for single band ID Mott insulators. The Hamiltonian is of the form 



H = —t 
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where n,- 



n j 4 are electron number operators and U > 0. 



2. ID Holstein modelS 3 -: 

The Holstein model provides an example of an incommensurate CDW state and describes a partially filled band 
of spin-1/2 electrons coupled to dispersionless phonons of frequency luq = 



H = -t 
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Integrating out the phonons induces a retarded attractive electron-electron interaction. In the limit t <C wo the 
retardation effects can be neglected, leading to an effective attractive Hubbard model with U cx —A 2 . 

3. Su-Schricffcr-Hceger model 2 ^: 

A second example of an incommensurate CDW state is provided by the Su-Schrieffer-Heeger model, which 
describes a partially filled band of spin-1/2 electrons coupled to dispersing phonons 



H = -Y i [t-KQj + i-Q i )] 
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Taking the continuum limit of ([5]) (which describes the behavior at low frequencies u < t) it was shown by 



Fradkin and Hirsch 2 ^ that the regime of large phonon frequencies t <C i^o 
Thirring mode&^. 



rr is described by the U(l) 



In all three cases a continuum description of the low-energy electronic degrees of freedom is obtained by considering 
only the modes in the vicinity of the Fermi points ±fcp. The lattice electron annihilation operators are expressed in 
terms of slowly varying right- and left-moving Fermi fields as 
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where do is the lattice spacing, x = jag, and a =t,i labels the spin. In the bulk the fields R a and L a are bosonizcd 
according to 

Rt(r,x) = 



L U T > X ) = ~j= ex p(-5^c(T,x)^) exp^-^/ CT ^ s (r,a;)), 



(7) 
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where the Klein factors r\ a satisfy anticommutation rules \r\a^f]'a\ = 2(5 CTfT ' and /-j- = 1 = — /j.. The fields a and </> a 
are the chiral components of the canonical Bose fields $ a and their dual fields 0„, 



$a = §a + 4>a, 6 a = 4>a - 4>a, O = C, S. 

In the bulk the Hamiltonian density then can be cast in the spin-charge separated form 
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The charge and spin velocities v C B , the Luttinger parameters K CjS and coupling constants g CjS are functions of the 
hopping integrals and interaction strengths defining the underlying microscopic model. The cases discussed above 
correspond to the following parameter regimes: 

1. Mott insulators 

As a result of repulsive electron-electron interactions we have v c > v s and K c < 1. The cos<f> c perturbation 
in the charge sector is relevant and opens up a gap. The cos $ s interaction in the spin sector is marginally 
irrelevant and flows to zero under the renormalization group. We therefore neglect it in the following. 

2. Electron-Phonon Systems 

At low energies the electron-phonon coupling induces an attractive electron-electron interaction, which results 
in v s > v c and K c > 1 > K s . The cos$ c term is irrelevant, while cos$ s term is relevant (marginally relevant in 
the spin-SU(2) symmetric case K s = 1) and opens up a gap in the spin sector. 

In both cases we end up with a spin-charge separated theory of a gapless Luttinger liquid and a sine-Gordon model. 

We now imagine a strong, local potential to be present. It is well known from the work of Kane and Fisher— that 
the coupling to the impurity is relevant, leading to a crossover at a characteristic dynamical energy (and temperature) 
scale called "Ik" , below which the system is effectively cut into two disconnected parts, and to a pinning of the CDW. 
This potential could be due to a strong potential impurity or an STM tip. We model the strong impurity potential 
by a boundary condition on the continuum electron field 



V a (x = 0) = Q. 



(11) 



An important physical consequence of the pinning of the CDW at the impurity (or boundary) is the development 
of an induced static CDW order in this (effectively) quantum critical system. This induced static CDW order, often 
referred to22 as a "Friedel oscillation" , leads to non-dispersive features in the LDOS^ which can be detected in STM 
and STS experiments. 

Our effective low-energy Hamiltonian in the CDW case (in the case of a Mott insulator the roles of spin and charge 
sectors are interchanged) then becomes 



H = H ^ 
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where the Bose fields are subject to the hard- wall boundary conditions (we consider more general boundary conditions 
in Sec. E} 



$cs(a; = 0) = 0. 



(15) 
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In Appendix [X] wc argue that a weak potential impurity rcnormalizes to strong coupling even for moderate attractive 
interactions, suggesting that this situation too can be modeled in terms of the boundary conditions (fl~5|) . We note 
that our starting point (fT2|) -(fl~5 j) differs from the model considered in Ref. [H, where the impurity couples only to the 
gapless charge sector. 

The charge sector (|13[) describes gapless collective charge excitations propagating with velocity v c , which carry 
charge =Fe and are commonly referred to as holons and antiholons respectively. On the other hand, the spin exci- 
tations or spinons are described by the sine-Gordon model on the half- line ([M]). which is known to be intcgrable 
for quite general boundary condition a 29 ' 30 . In the regime K s > l/y/2 the elementary bulk excitations are gapped 
solitons and antisolitons which correspond to up- and down-spin spinons respectively. For K s < 1/ s/2 propagating 
breather (soliton-antisoliton) bound states occur as well. At the Luther- Emery point (LEP) K s = l/y/2 the spin 
sector is equivalent to a free massive Dirac fermion 3 -^. The exact bulk scattering matrix was first derived by Zamolod- 
chikov^; the boundary reflection matrices of solitons and antisolitons 2 ^ and breathers^ 3 , were derived by Ghoshal and 
Zamolodchikov. We will restrict ourselves to the regime A^ s > l/y/2 throughout, which implies that no breathers 
exist. 

The lattice models discussed above give rise to the simplest kind of one dimensional CDW state/Mott insulator. 
More complicated versions arise in strongly correlated two- and three-leg ladder system s 19 ' 34 . In these systems, even 
though electron interactions are strongly repulsive, a Mott state with a finite (and typically large) spin gap is found 
for a range of dopings close to half filling. While the precise description of the spin sector for two-leg ladders is 
considerably more complicated 3 ^, we expect our calculation to capture important qualitative features. 



III. GREEN FUNCTION 



The central object of our study is the time-ordered Green function in Euclidean space, 



G aa '(T,X 1 ,X 2 ) 
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where | Ob) is the ground state of (|12|) in the presence of the boundary and r = it denotes imaginary time. The spin 
takes the values f and J,. At low energies the linearization around the Fermi points yields the decomposition 
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where e.g. G R ^, = — (0b| T T R(t(t, Xi) ^v(0, 2:2) |0b). As we are interested in the LDOS, we ultimately want to set 
X\ = X2- Below we will calculate the spatial Fourier transform of the LDOS as physical properties can be more easily 
identified. In momentum space the RL and LR contributions occur in a different region (Q ss ±2fcp) compared to the 
RR and LL parts (Q ~ 0). In absence of a boundary we have G R ^, = G^ R = as the charge parts of these Green 
functions vanish. In presence of a boundary left and right sectors are coupled and the Fourier transform of the Green 
function (|17[) concomitantly acquires a nonzero component at Q w ±2fcp, which provides a particularly clean way of 
investigating boundary effects. For this reason we first focus on the 2/cp-part of the Green function and study the 
small momentum regime afterwards. 

The Green function G R ^, factorizes into a product of correlation functions in the spin and charge sectors. The 
charge part can be determined by standard methods 3 ^— (see App. [B]). On the other hand, the intcgrability of the 
sine-Gordon model on the half-line (|14p enables us to calculate correlation functions in the spin sector using the 
boundary state formalism introduced by Ghoshal and Zamolodchikov 2 ^ together with a form-factor expansion 3 ^—. 
As we show in App. [C]t,he leading terms of this expansion yield (r > 0, Xi < X2) 



g c (T,xi,x 2 ) 



i,x 2 ) = g c (r,xi, 

bar,' 



X2)9b(t,Xi,X 2 ), 

1 



1 



2- 



2iR) 



2iR) 



4x]_x 2 



(v c t — ir)(v c T + ir) 



g s (T,xi,x 2 ) = Zie 1 * 



■tt\ „9/2 2i— flsinhS -At cosh , 
1 2 J c fc- s fc- "1 r 



(18) 
(19) 

(20) 



where g CtS arc the contributions of the charge and spin sectors respectively. Here Kq is a modified Bcsscl function and 
the center-of-mass coordinates are R = (x\ + x 2 )/2 < and r — x\ — x 2 < 0. The normalization constant Z\ was 
obtained in Ref. [42|. At the LEP— and the SU(2) invariant point the boundary reflection amplitude K{6) is given 
by 



K(0)=i tanh - for K s 
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The expressions for general values of K s can be found in Rcfs. 
related to the Luttinger parameter by 



[29ll46ll47j . The exponents in the charge sector arc 
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We stress that (J2DJ) is independent of a and that the dependence on K s is through the overall normalization constant 
Zi and the boundary reflection amplitude only. The one-particle contributions of the form-factor expansion are given 
by the first two terms in (|20[). while the dots represent corrections involving a higher number of particles in the 
intermediate state as well as higher-order corrections due to the boundary. We have determined the subleading terms 
in the spin part of the Green function at the LEP and found their contribution to the LDOS calculated below to be 
negligible (see Sec. lIVC]) . 

After analytic continuation to real times r — >• it the Green function (fl8]l exhibits a light cone effect. The first term 
in (|20p shows oscillating behavior for r 2 < (v s t) 2 but is damped otherwise. Similarly, the second term is oscillating 
for 4i? 2 < (v s t) 2 . These oscillations are due to the propagation of spinons from x 2 to x\ either directly or via the 
boundary. In particular, at late enough times both terms will oscillate. A similar light cone effect was observed^ in 
the Ising model with a boundary. On the other hand, the charge part (|19p possesses singularities at v c t — ±r and 
v c t = ±2R due to the propagation of antiholons. 

The small momentum regime of the Fourier transform of the LDOS is obtained from (r > 0, x\ < x 2 ) 
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Compared to G RL the singularity at v c t — 2R is much softer whereas the one at v c t — r is more pronounced. 



(23) 



IV. LOCAL DENSITY OF STATES 

The knowledge of the Green function (fTSf enables us to calculate the LDOS, which is directly related to the tunneling 
current measured in STM experiments. As was noted in Ref. it is useful to consider the Fourier transform of the 
LDOS as physical properties can be more easily identified. For example, this technique was used to study quasiparticle 
interference in high-temperature superconductors^ and rare-earth compounds^. We will first consider the boundary 
condition ^(O) = 0, which results in the Green function (JT5J) and hence yields a spin-independent LDOS. More 
general boundary conditions may lead to a spin-dependent LDOS or even the formation of a boundary bound state. 
We will discuss this case the next section. 

The Fourier transform of the LDOS is given by N a (E, Q) = N>(E, Q) + N<(E, Q), where 



N>(E,Q) = - — 



N<{E,Q) = — 
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Here the Green function has been analytically continued to real times and we have take the limit x\ — > x 2 = x. We 
will focus on the LDOS for positive energies in what follows but note that the LDOS for negative energies can be 
analyzed analogously. As mentioned before, we will be mainly concerned with the 2/cF-component as it vanishes in 
the absence of the boundary and hence offers a particularly clean way of investigating boundary effects. For Q « 2kp 
only contributes and starting from (|T8|) we arrive at the following expression (see App. [D| 
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FIG. 1: \N& (E,2kF + q)\ (arbitrary units) for K c = 0.8, K s = 1 and v c = 2v a . The curves are constant g-scans which have 
been offset along the y-axis by a constant with respect to one another. The LDOS i s dominated by a strong peak at q = 0, 
i.e. Q — 2&f. We further observe dispersing features at E c = v c \q\/2 + A and E s = \J (v s q/2) 2 + A 2 . For q < the dispersing 
features are strongly suppressed. 

Here \q\ <C 2kp, A = arcosh(^), F\ denotes Appell's hypergeometric function^ (see App. [E]), h\(ff) = 1, Ii2(9) = 
A'(6» + if) e 9 / 2 , and 

2 / . , „\ . ^v s q\ r 2v s E — Acosh.6 . fv s q „ . , a . 

«i = — ^- Acosh^ +isgn (-p U, u 2 = ^ oa ■ , fl +isgn(-^-2sinhg)j, 28 

v c q V A / v c V s q — 2Asmn0 V A / 

where (5 -> 0+. The result ((27j) is valid for a + & < 2 and -1/2 < c. Below we plot N>(E, 2&f + g) for two different 
parameter regimes. We smear out singularities by taking 5 small but finite (6 = 0.01 unless stated), which mimics 
broadening by instrumental resolution and temperature in experiments. The results presented below apply to the 
regime T <C E, A, v c /ao (arj is the lattice spacing), where temperature effects are negligible. 

A. Repulsive Case 

We first consider the case v s < v c , K c < 1. This can be thought of as providing a simplified model for the LDOS of 
a two-leg ladder with repulsive electron-electron interactions. The low-energy theory for the ladder is similar in that 
there is a gapless charge sector and a gapped spin sector, but the full description of the latter is considerably more 
complicated^. 

In Figs. [1] and [2] we plot N>(E, 2kp + q) for the case of unbroken spin rotational symmetry {K s = 1). The Fourier 
transform of the LDOS is dominated by a singularity at momentum 2fcp (q = 0), which arises from the contribution 
Ni . For fixed energy and close to the singularity this term behaves as (see App. |F]) 

N>(E,2k F + q)~ ( — ) , a=l-^, (29) 
\v c qj 2 

which implies a phase jump of ira as q — ¥ 0±. This peak is indicative of the CDW order being pinned at the boundary. 
We note, however, that the peak occurs at finite energies and hence the underlying process is not static. A similar 
feature can be seen in the Luttinger liquid case£, where the singularity as a function of q is softer (a^L = (1 — K 2 )/2). 

At low energies above the spin gap A we further observe two dispersing features associated with the collective spin 
and charge degrees of freedom respectively. These are broadly similar to the bulk single-particle spectral functio n 50 ^ 1 
and feature (1) a "charge peak" that follows 

E c (q)= V -^ + A (30) 

and (2) a "spin peak" at position 
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FIG. 2: Constant energy scan for E = 2A: \N>(E,2k F + q)\ (arbitrary units) and ArgN>(E,2k F + q) for K c = 0.8, K „ = 1 
and v c = 2v s . We observe a peak at q = (related to the pinning of the CDW at the boundary) and dispersing features at 
q = ±2(E — A)/v c as well as q = ±2V E 2 — A 2 /v s . For g < the dispersing features are strongly suppressed. Furthermore, we 
observe characteristic jumps in the argument at the positions of the peaks. 



We note that neither peak is sharp (i.e. they are not delta-functions) and hence have to be thought of as arising from 
excitations involving at least two "elementary" constituents. In this way of thinking, the charge peak arises from 
two-particle excitations composed of a "zero momentum" spinon contributing an energy A and a gapless antiholon of 
"momentum" q. On the other hand, the spin peak can be thought of arising from two-particle excitations composed 
of a "zero momentum" antiholon and a spinon of "momentum" q. The appearance of v c /2 and v s /2 in (|30[) and 
(|31[) . respectively, is due to the fact that the particles have to propagate to the boundary and back, thus covering 
the distance 2x in time t. We note that on a technical level the charge peak arises from the contribution N{~ to the 
Fourier transform of the LDOS, whereas the spin peak has its origin in N^", which encodes the effects of the boundary 
on the spin degrees of freedom. In the q < region the dispersing features are strongly suppressed and for K c = 1 
the charge feature is found to vanish entirely. 

It is instructive to plot N>(E, 2&f + q) as a function of q for fixed energy, see Fig. O We observe characteristic 
jumps in the phase Arg at the peak positions. This is similar to the Luttinger liquid case^. 



B. Attractive case 



We now turn to the case of a CDW state arising in a system with (effective) attractive electron-electron interactions. 
As discussed in Sec. [HI this case arises in electron-phonon systems. The effective parameters are given by v s > v c and 
K c > 1 > K s . 

In Fig. [3] we plot N> (E, 2kp + q) as a function of energy for several values of q (in units of the spin gap). We again 
observe a singularity at 2fcp, which arises from (|29p . The singularity is much less pronounced than in the repulsive 
case and disappears for K c > y/2. Like in the repulsive case there are several dispersing features: 

1 . a charge peak at E = E c (q) , where E c is given by ([30]) ; 

2. a spin peak at E = E s (q), where E s is defined in (pTTj) ; 

3. when \q\ exceeds a critical value go a third dispersing low-energy peak appears (see Fig. 2]) at 



E cs (q) = V -^- + Adl- 



\ll-( V -7-) 2 =EM + V f(\ q \-qo), qo = —^=. (32) 



This feature can be thought of as arising from a "momentum" qo spinon and an antiholon carrying "momentum" 
q — qo- We note that in this case the spin and charge excitations have the same group velocity 

dE c dE s v c , . 

dq dq q=qo 2 



This behavior is reminiscent of what is found for the single-particle spectral function in the bul k 50 ' 51 . The peak 
splitting and hence the qualitative difference between the repulsive and attractive regime is a consequence of the 
curvature of the (anti)soliton dispersion relation and hence of the spin gap. In the Luttinger liquid case^ (where both 
sectors are massless) there are only two dispersing features in both regimes. 
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FIG. 3: |JV^"(S, 2/cf + q)| (arbitrary units) for K c = 1.2, if s = 1 and v s — 2v c . The curves are constant g-scans which have 
been offset along the y-axis by a constant with respect to one another. The peak at q — is much less pro nounced then in the 
repulsive case (see Fig. [T]). We observe dispersing features at E c , E s , and E cs = v c \q\/2 + Ayl — (v c /v a ) 2 (for \q\ > go only). 




FIG. 4: \N>(E,2k F +q)\ (arbitrary units) for K c = 1, K a = 1/^2, v B = 1.5 v c , S = 0.001, and v c q/A = 1.0, . . . , 1.8. The curves 
have been offset along the y-axis by a constant with respect to one another. We observe the splitting of the spin peak at E a at 
the critical momentum v c qo/A ss 1.19. 



C. Higher-order corrections 



As we have indicated in Eq. ([20)1 there are contributions to the LDOS beyond those that we have discussed above. 
They arise from our calculation of the spin-part of the Green function and are expected to be small^. In order to 
verify that they can indeed be neglected, we have analyzed them in some detail at the LEP K s = 1/V2, where the 
necessary matrix elements take a particularly simple form, which makes the actual calculations much easier. 

Our purpose is then to determine further terms in the the expansion l[26|) of N>(E 7 2kp + q). We denote by N nm the 
contribution to (|26j) that arises from processes in which n gapped spinons (which correspond to solitons/antisolitons 
in the sine-Gordon model describing the spin sector) propagate between (0,cc2) and (r, x\), and which involves the 
m'th power of the boundary reflection matrix K. The terms discussed above correspond to TVio = N{" and Nqi = TV^ . 
In App. [C4l we present some details on the calculation of the terms in the form-factor expansion of <?s(t, xi, X2) that 
give rise to N nm with m + n < 3. The number of ^-integrations in N nm equals m + n (cf. (|27j) ). We find that their 
contributions to the Fourier transform of the LDOS are small. In particular, all qualitative features of the LDOS 
such as dispersing peaks are already encoded in N\q and Nqi . In Fig. [5] we show the leading terms N\q and Nqi as 
well as the sub-leading terms for K c = 1 and v c > v s for a fixed value of q as a function of energy. We see that the 
two leading terms in (|26|) indeed capture all qualitative features of the LDOS and carry the main part of the spectral 
weight at low energies E < 5 A. The higher-order terms are small compared to N\q and Nq\. In particular, N30 
vanishes for E < 3A, since this term originates from a three-particle process. In general, all terms N nm originating 
from n-particle processes vanish for E < nA. Most importantly, however, the higher-order terms do not possess 
any singularities. The suppression of subleading terms in the form-factor expansion for bulk two-point functions is a 
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FIG. 5: (Color online) Left: Comparison of the absolute values of Nw — N{"(E,2kY + q) and Noi = A^(-E, 2&f + q) to the 
sub-leading terms Nu and N02 (N20 = 0). We stress that the scale for the higher-order terms has been magnified. Right: 
Absolute values of the terms N30, N21, N12, and N03. We stress the different scales on the y-axis. The parameters are K c = 1, 
K B = 1/V2, v c — 2v B , and v B q/A = 3. The three-particle contributions N30 and N21 vanish for E < 3A. Furthermore, the 
higher-order terms possess no peaks. 



well-known feature of massive theorie o ' 52 i 53 , whereas the smallness of terms involving higher powers of the boundary 
reflection amplitude K has recently been demonstrated for the Ising model with a boundary magnetic field^. 



D. Small-momentum regime 



The small-momentum regime Q as of the Fourier transform of the LDOS can be analyzed in the same way as 
in the Q ps 2/cp case discussed above. We note that the LDOS for Q as is non-vanishing even in absence of a 
boundar y 50 ! 51 . In the presence of a boundary the Fourier transform of the LDOS for Q as is obtained from (|23| . 
The leading terms are given by N>(E, Q) = Xa=i N?(E, Q) + . . ., where 

Z, p-if ( 2c +!) r( a + h+l) 
N>(E,Q) = -B(E-A) Zl 



4tt 2 v2 +b - 1 r(2a + 2b) 

M^K 2c+1 



A V-Aco S hV-^ F^ + l,c,c,2a + 26; <, (34) 

Here we have \Q\ <€. kp, A = arcosh(^), h\{6) = e 9 / 2 , h,2(0) = K(6 + V^\, and u\ t 2 are defined in with q replaced 
by Q. The main difference to (|26|) is in the dependence of Appcll's hypcrgcometric function on K c . 

In Figs. |5]and [7] we plot N>(E,Q) for the case of repulsive electron interactions and unbroken spin rotational 
symmetry. It is dominated by a singularity at Q = 0, which has its origin in N{" and behaves as ~ l/Q independently 
of K c . This singularity is more pronounced than its counterpart at 2fcp- We further observe dispersing features at 
positions E C (Q) and E S (Q) respectively. Both of these are symmetric under Q —> —Q. The peak at E C {Q) is strongly 
suppressed, vanishes for K c = 1, and becomes a dip in the attractive regime. The suppression is due to the softness 
of the singularities of G RR at v c t = 2R. On the other hand, the charge part of (|2"3"]) has its strongest singularity at 
v c t = r = 0, which results in a background of spectral weight in N> (E, Q) for all energies above the spin gap. 

In the attractive case (v s > v c ) we observe a peak-splitting similiar to that in the 2fcp-component, but all peaks are 
very weak. 



V. GENERAL BOUNDARY CONDITIONS AND BOUNDARY BOUND STATES 



So far we have considered the simplest possible boundary conditions corresponding to a spin-independent phase 
shift of 7r. Both ways of realizing a boundary in a (quasi) one dimension system that we have discussed above (i.e. as 
a result of an impurity or in a "two-tip" STS experiment) are expected to give rise to a local potential or magnetic 
field. These correspond to more general phase shifts for reflection of particles at the boundary. As is well known, 
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FIG. 6: \Na(E,Q)\ (arbitrary units) for K c — 0.8, K s — 1 and v c = 2v s . The curves are constant Q-scans which have been 
offset along the y-axis by a constant with respect to one another. We observe dispersing features at E C (Q) and E S (Q). The 
charge features are very weak for all momenta. We further observe a constant background at energies larger than the spin gap. 




FIG. 7: Constant energy scan for E = 2A: \N>{E,Q)\ (arbitrary units) and ArgN>(E,Q) for K c = 0. 8, K s = 1 and v c = 2v B . 
We observe a strong peak at Q = and dispersing features at Q = ±2(E — A)/v c as well as Q = ±2\/ E 2 — A 2 /v s . Note that 
\Na (E,Q)\ is symmetric under Q — > —Q. 



such more general boundary conditions can give rise to boundary bound states, see for example Refs. [54j . These 
are expected to be visible in the Fourier transform of the LDOS as "resonances" inside the single-particle gap. This 
is most easily seen by considering a Lehmann representation of N>(E,Q) in terms of the cigenstates |rib) on the 
half-line 

r° 

N>(E,Q)= dxe- iQx J2 \(0 h \^a(0,x)\n h )\ 2 S(E-E nb ). (35) 
J -°° 

For boundary bound states |bbs, a) we have < E'bbs.a < A, which leads to features in N>(E,Q) below the single- 
particle gap. As we are dealing with a spin-charge separated system, these features will generally not be sharp as the 
bound state occurs only in the gapped sector of the theory. We now turn to calculating the LDOS in cases where 
boundary bound states exist. We first consider boundary conditions of the form 

R a (r,0) = -e"^*°/ 2 Mt,0), (36) 

where /-)- = 1 = — f^. In terms of the Bose fields these boundary conditions read 

$ c (r,0)=0, $ s (r,0) = $ s °, 0<$ S <7T. (37) 

We note that these boundary conditions break spin rotational symmetry. However, if we go over to the case of a 
one-dimensional Mott insulator by exchanging spin and charge degrees of freedom, the spin rotational symmetry 
remains intact and the boundary conditions correspond to a local potential. 

As before we will focus on the 2fcp-component of the Fourier transform of the LDOS. As wc have changed the 
boundary conditions only in the spin sector, the charge part (|19[) of the chiral Green function remains unchanged. 
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The two leading terms of the form- factor expansion in the spin sector are still of the form ([20|) . but now the 
boundary reflection amplitude K is different and in particular is spin-dependent. At the LEP it is given by^ 5 - 

sin (if - /<r%) 

K° s {6) = — )- j^-, K aa {9) = K ss {9) = 0. (38) 

C09^i 3 + /«r-f J 

Here we have introduced the notations ~\— + and 1= — as well as a = — for a = + and vice versa. We note that as a 
result of a different choice of phase for the asymptotic states E q. (l38l) differs by a minus sign from Ref. [45| (see also 
App. [Cj) . The expressions for general K s can be found in Refs. |29ll4al47j . 

In the spin rotationally symmetric case K s — 1 we have (9) = K ^(ff) = K(0), where K(9) is given in (|21l) . 
We stress that the Green function remains diagonal in spin space, G^, oc S aa >, and that the spin-dependence is 
entirely due to the boundary reflection matrix K aa . Before presenting the resulting LDOS we discuss the emergence 
of a boundary bound state in the spin secto r 29 ' 46 ! 55 . If we choose the phase-shift in the spin sector such that 

Kl-K < <&°, (39) 

the boundary reflection amplitude K h (#) has a pole in the physical strip < 3m 9 < ir/2. This pole corresponds to 
a boundary bound state with energy 

£ bbs =Asin 7 , 7 = (40) 

The physical nature of the bound state has been discussed by Ghoshal and Zamolodchiko\«2 9 .. The classical ground 
state of the sine-Gordon model on the half-line (fT4"f in the entire range < <&g < tt is characterized by the asymptotic 
behavior $ s — > as x — > — oo. On the other hand, there exists a second classically stable state satisfying $ s — > 2tt as 
x —> — oo. When is sufficiently large this state is expected to be stable in the quantum theory as well. 

We note that for $!? = it both states are degenerate and vanishes. In the attractive regime of the sine-Gordon 
model K s < l/y/2 additional boundary bound states occur, while in the spin rotationally invariant case K s — > 1 the 
condition it < < ir is never satisfied and hence no boundary bound states exist. 

When calculating dynamical response functions in the boundary state formalism additional contributions in the 
form factor expansions occur upon analytical continuation in the rapidity variables. In particular, the pole of the 
boundary reflection amplitude in the physical strip gives rise to an additional term linear in K in the form-factor 
expansion ([2"U]) . In the case r > and x\ < x-i it takes the form (see Add. IC 5"|) 

e($°-ir s 2 7r) 5 al Z 1 Be^ e 2 ^ Rcosi e- Arsin \ (41) 

where the constant B > is related to the residue of K T± (9) (see f[C47J> ). At the LEP it equals B = -2cos$|?. We 
stress that this additional term appears in the down-spin channel only, since we have assumed < $g < tt. If we were 
to consider — tt < $^ < 0, we would find a term similar to (|41 [) in the up-spin channel only. 

The Fourier transform of the LDOS for the boundary conditions (|37p can be expanded as before and is expressed 



N>(E, 2k F + q)=Y, N>4 E ' 2fcp + 9 ) + . . . . (42) 

Here the first two terms are again of the form (f27|) . where in the second term N^ 2 we need to replace the boundary 
reflection amplitude K(d + \^) by its spin-dependent counterpart K aa {9 + \^). The third term is obtained from (j4T|) 
and arises as a result of the presence of a boundary bound state. Explicitly it reads 



2c+l 

(3) 



(E — -Ebbs) 



F^> (2c + 1, a, 6, 2c, a + b + 2c; «*, -u 3 , -u' 3 ) . (43) 



Here denotes Lauricella's hypergeometric function of three arguments^ 6 - (see App. [E| and 

2 2A 2A 

u 3 = (E-E hhs )+i — COS7, u' 3 = i COS7, (44) 

v c q v s q v s q 
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FIG. 8: \N^ (E,2kp + q)\ (arbitrary units) for K c = 1, K s = 1/V2, v c = 2v B , and $° = 0.97T. The curves are constant q-scans 
which have been offset along the y-axis by a constant with respect to one another. We observe dispersing features at E s (q) and 
E c (q) (for q > only) as well as a non-dispersing singularity at E = -Ebbs, which is due to the formation of a boundary bound 
state in the spin sector. 




E/A 

FIG. 9: \N>{E,2k F + q)\ (full line) and \N>(E,2k F + q)\ (dashed line) for v s q/A = 6, K c = 1, K a = 1/V2, v s = 2v c , and 
$° = 0.9 7T. The broad maximum at .E ~ 1.8 is caused by the excitation of the boundary bound state and additional charge 
excitations. 



where the constant 7 defined in (|4"U)) . The Fourier transform of the LDOS ((32) has a non-dispersing singularity at its 
lower threshold 

NUE^ + ^cc^-^, a=l-^. (45) 

The emergence of a non-dispersing feature within the spin gap signals the presence of a boundary bound state. In 
Fourier space the LDOS is a convolution of contributions from the spin and charge sectors. As we are dealing with a 
bound state in the spin sector, the exponent of the singularity depends only on the Luttinger parameter in the charge 
sector. We note that the singularity occurs only in the down-spin channel and disappears for < 1/2. On the other 
hand, in N< the additional feature due to the boundary bound state appears only in the up-spin channel. 

In Fig. [8] we plot the down-spin component of (|42|) for v c > v B as a function of energy for several values of q. As 
before, at low energies above the spin gap A we observe two dispersing features associated with the collective spin and 
charge degrees of freedom that follow E s and E c respectively. In addition, we observe the non-dispersing singularity 
gl]) at E = E hhs . 

In Fig. [S]we plot N> and as functions of energy for v c < v s . We see that the singularity arising due to the 
presence of a boundary bound state appears only in the down-spin channel. For either spin polarization we observe 
three dispersing features at E c (q) 7 E s (q) and E cs (q) respectively. Their interpretations are completely analogous to 
discussion in Sec. IIV Bl In addition to these sharp peaks we observe a broad maximum in the down-spin channel at 
energies E w Ebbs + v c q/2. This feature is suppressed for v c > v s , see Fig. [H Its physical origin is the simultaneous 
excitation of the boundary bound state and a finite energy excitation in the charge sector. We note that the asymmetry 
in (E, Q) — N> (E 7 Q) could in principle be detected in experiments using a magnetic STM tip. So far we have 
considered only hard- wall boundary conditions in the charge sector, i.e. $> c (x — 0) = 0. Our analysis can be 
straightforwardly extended to the case $ c (^ = 0) = 3>°, which in terms of the original electrons corresponds to a local 
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FIG. 10: \N>(E, 2k F + q)\ (arbitrary units) for v s q/A = 2, K c = K s = 1, and v c = 2v a . We observe spectral weight within the 
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FIG. 11: Constant energy scan of \N^ (E, 2/cf + q)\ (arbitrary units) for E — 2A, K c = K s = 1, and v c = 2v s . We observe a 
suppression of the peak at q = related to the pinned CDW. 



potential close to the boundary. 



VI. FINITE-TEMPERATURE LDOS 



Another interesting issue concerns the effects of a finite temperature on the LDOS. The regime T < A/2 can in 
principle be analyzed by generalizing the methods recently developed in Refs. [57| to the boundary state formalism. 
However, in order to keep matters simple we will restrict ourselves to the regime of very low temperatures T <C A. 
Here the main effects arise from a modification of the dynamical response in the gapless charge sector and correlation 
functions in the spin sector can be approximated by their T = expressions. This is the case because we only 
consider response functions that involve both sectors. The charge part of the Green function G^,(t,x\,X2) = 
5 c (r, x\,x 2 ) 5s(t, xi, x 2 ) can be evaluated using conformal field theory method o 38 ' 58 and is found to be 



g c (r,xi,x 2 ) 



U (J<J' 



7T 



a+b 



sm a (^(v c r - 2iR)) sm b {^(v c r + 2LR)) 

off' 



sin (ffl( w c^- irj) sin(^(v c r + ir)) 



(46) 



Here (3 = l/k^T and the exponents in the charge sector were defined in (|2"2"1) . As we have already discussed, the spin 
part g s (r,xi,x 2 ) is given by ([20]). 

The particle contribution to the Fourier transform of the LDOS is still given by ([24]) . The form- factor expansion 
in the spin sector results in the series expansion N>(E, 2kp + q) = J2i N(*(E, 2kp + q), where the first two terms can 
be cast in the form 



N>(E,2k F + q) 



ZiTT a+b e'T e _i ^ ( a + b + 2c ) 



32tt 5 w c (v c f3) a+b - 2 

poo 

d9hi(6) / dx 



dy 



( sinh(i(y - x)) 



sinh a (a: - id) smh b (y-i6) \ sinh(i(:r + y - id)) 



(47) 
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FIG. 12: \N> (E,2kp + q)\ (arbitrary units) for v B q/A — 6, K c = K s = 1, and v B — 2v c . We observe again that the broadening 
of the propagating spin peak is stronger than that of the other peaks. 



hi(0) = l, h 2 (6) = K a9 (8 + ii) e e/2 , E 1 = E 2 = E - A cosh 6, qi = q, q 2 = q sinhfl. (48) 

We have plotted N> (E, 2kp + q) for different temperatures in Figs. fTUHT^I As expected a finite temperature leads 
to a softening of the spectral gap A, a suppression of the peak related to the pinned CDW, and a broadening of the 
dispersing peaks. We observe that the effect of an increasing temperature on the spin peak is much stronger than the 
effect on the charge peak. The physical reason for this is as follows: In the CDW state only the spin sector is protected 
by the gap. Thus for T <C A there exists a significant number of antiholons in the thermal ground state. They will 
participate in the distribution of the external momentum q after the creation of an additional antiholon-spinon pair, 
thus leading to a decreased probability for the spinon to take the momentum q and thus to a suppression of the spin 
peak following E s (q). On the other hand, the charge peak is not affected as the antiholons possess a linear dispersion. 
This behavior is reminiscent of what is found for the bulk spectral functions^. 



VII. IMPLICATIONS FOR STM EXPERIMENTS 



STM experiments measure the local tunneling current, which is related to the LDOS by Eq. ([T]). In particular, 
the voltage dependence of the tunneling conductance measures the thermally smeared LDOS of the sample at the 
position of the tip. A possible spin dependence in the LDOS can be detected using a magnetic tip. As we have 
considered a one-dimensional model, our results apply to quasi-lD materials at energies above the 1D-3D cross-over 
scale, which is set by the strength of the 3D couplings. Furthermore, the main feature of the model we have studied 
is the existence of a spectral gap in one of the sectors, while the excitations in the other sector remain gapless. This 
situation is experimentally realized in various materials, for example in two-leg ladder materials^, stripe phases of 
FITSC^, carbon nanotubesi^, Bechgaard salts^, and chain materials^ like SrCu02 and S^CuOa. As our results 
show, STM experiments can be used to extract rather detailed information regarding bulk excitations by analyzing 
the modification of the LDOS due to a boundary/impurity. 

Perhaps the most interesting materials to which our findings may be applied on a qualitative level are two-leg 
ladders liko^ Sri4Cu2404i which possesses a spin gap££ of Aw 550K. The model we have studied captures the most 
basic features of the low-energy description of (weakly doped) two-leg ladders, namely a gapless charge sector and a 
gapped spin sector. While the description of the spin sector for weakly doped two-leg ladders is considerably more 
involved, we expect the gross features to be similar. In particular, we expect peaks to appear in N a (E,Q), which 
correspond to the pinned CDW order, dispersing spin and charge degrees of freedom, and possibly boundary bound 
states. We note that our results apply to the regime T<£,A, v c /ao, where temperature effects are negligible. 

In Ref. @ it was proposed that STM and STS experiments in HTSC can be used to detect "fluctuating stripes" 
(i.e. incommensurate spin and charge fluctuations on energy scales small compared to the superconducting gap) by 
rendering them static by the effects of impurities with a potential comparable to the (low) energy scales of these 
fluctuations. In that work it was also argued that ID Luttinger liquids are effectively quantum critical systems 
and that a form of local (power law) CDW order is effectively induced by impurities (and edges) which pin the 
phase of the CDW, an effect that we have shown here to take place in a ID Luther-Emery liquid associated with a 
CDW. STM and STS experiments in E^S^CaC^Os+a have confirmed the existence of both non-dispersive spectral 
features in the LDOS associated with "fluctuating stripe order" as well as dispersive features associated with the 
propagating quasiparticlcs of the superconductor^^. Recent STS experiments in E^Si^CaC^Os+a have shown that 
the dispersive features of the LDOS disappear above the T c of the superconductor while the non-dispersive features 
survive up to the temperature T* at which the pscudogap closes^. 
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VIII. CONCLUSIONS 



In this work we have determined the spatial Fourier transform of the LDOS of one-dimensional CDW states and 
Mott insulators in presence of a boundary. The latter may either model a strong potential impurity or be realized 
in a two-tip STM experiment. We found that the Fourier transform of the LDOS is dominated by a singularity at 
an energy equal to the single-particle gap A and at momentum 2fcp. This feature is indicative of the pinning of the 
CDW order at the position of the impurity. We observed clear signatures of dispersing spin and charge excitations, 
which can be used to infer the nature of the underlying electron-electron interactions. In the case of CDW states 
with repulsive interactions we find a spin mode and a linear dispersing charge mode while for attractive interactions 
a third dispersing mode appears, which can be thought of as arising from a spin excitation with a fixed momentum 
qo and a charge excitation with momentum q — q$. 

We have also investigated the modification of the LDOS due to boundary bound states. These may arise in 
presence of boundary potentials or magnetic fields. We found that boundary bound states give rise to non-dispersing 
singularities at energies below the single-particle gap. While the bound state is formed in the gapped sector of the 
theory, the exponent of the corresponding singularity only depends on the Luttingcr liquid parameter of the gaplcss 
sector. We have analyzed temperature effects in regime T <C A and discussed implications of our results for STM 
measurements on quasi-lD materials such as doped two-leg ladders. 
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Appendix A: Renormalization group analysis of an impurity potential 



Let us consider the low-energy theory of a one-dimensional CDW state on the infinite line H = H c 



H r . = 



H s = 



Vc 

16tt 

Va 

16tt 



dx 



dx 



1 



( dx <s> s ) 2 + K*(d x e s ) 



(27T) 



dx cos <f> s 



We want to study the effect of an impurity potential at position x = 0, which in bosonized form reads 



Mmp — A 



dxS(x) cos 



cos 



H s , where 

(Al) 

(A2) 
(A3) 



As the spin sector in the bulk is massive, we have (cos($ s /2)) ^ 0, which implies that at low energies we can 
approximate cos(<I> c /2) cos($ s /2) in (|A3jl by (cos($ c /2)) cos($ s /2) + cos($ c /2)(cos($ s /2)). Thus in the charge 
sector we get a boundary sine-Gordon model^. For K 2 < 2 the impurity scattering potential scales to strong 
coupling. Hence as long as the interactions are not too attractive the field $ c gets pinned at the boundary, $ c (0) = 0. 
This in turn induces an impurity contribution in the gapped spin sector 

r dxS(x) cos 



V imp>s = A < COS 



(A4) 



If one analyzes the bulk and boundary cosine terms in the resulting impurity model (|A2[) and (|A4[) simultaneously, 
the leading order renormalization group equations are given by the scaling dimensions of the perturbing operators 



x. 



(A5) 



As long as K 2 > 2/3 the boundary term grows more rapidly than the bulk term. Assuming that it reaches the 
strong-coupling regime first leads to the pinning of the spin field "^(O) — 0- This cuts the chain in two half-lines and 
we obtain the model ([X2 ^) -([I5 ]) . 
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Appendix B: Calculation of the Green function: Charge sector 



The Green function (JT6J) factorizes into a product of correlation functions in the spin and charge sectors. For 
example, using the bosonization identities ([7|) and ((HJ) the 2fcp-component G^, can be written as 

Both correlation functions have to be determined in the presence of the boundary at x = 0. The charge part is 
calculated below using a standard mode expansion^, the spin part will be calculated in App. [Cj 

In order to obtain the correlation functions in the charge sector we first bring the Hamiltonian (|13[) to standard 
form by rescaling the fields as $ c — > A' c $ c , C — > 9 C /A' C . The charge parts of the operators (JT)) and then become 



(±^ c (r,x)) -> e i " sc/4 exp(±ic0 c (z))exp(±i S c (z)) =e- i ^ c/4 exp(±i S c (z))exp(±i C( /) c (z)), (B2) 
(±i0 c (r,x)) -> e iffsc / 4 exp(±is^ c (z))exp(±i C c (z)) = e " i7rsc / 4 cxp(±i C ^ c (z)) cxp(±i S c (z)) , (B3) 



exp 
exp 

where we have already used (|B7[) and assumed — L < £ < 0. The constants are parameterized via s = sinh£ c and 
c = cosh£ c with K c = e". The complex coordinates are defined as z = v c t — ix,z = v c t + ix. The charge part of the 
Green function can hence be obtained from the four-point function 

gij8i0 o (*i) giai0o(zi) e ia 2 0c(z2) e i/3 2 c (2 2 )\ (B4) 

I UHP' 

where 0:1,2, 2 £ 1, zi = ij c t — vx\ and Z2 = — ia; 2 lie in the upper half-plane. 

We calculate (|B4I) from the mode expansions for the chiral fields cj) c and </> c . These are obtained by first noting that 
the fields $ c and C have to satisfy the equations of motion v c d x <d c = — i<9 r $ c and d T Q c = iv c d x & c as well as the 
boundary conditions $> c (% = 0) = & c (x = —L) = 0. The semi-infinite system is obtained by taking L — > oo. This 
yields the mode expansions 

OO • TL7TX 

$c(t,x) = -fio+iV (b n e- n ™^' L - bl e n ™°T' L ) , (B5) 

n—1 v 

OO TL7ZX 

e c (r, x) = 9o - i ^n + V (&„ e-™«W £ + fet e -cr/xj (B6) 



where the zero-mode operator IIo has the discrete spectrum 2nm, m £ Z, and 



e ,n 



87ri, r& m ,6^1 = 87T(5„ 



The mode expansions for the chiral fields are easily obtained via <p c = ($ c + O c )/2 and </> c = ($ c — c )/2. Their 
commutation relations are 

, a; = x' = 0, 

[<t> c (T,x),<l> c (T,x')] = { 47ri , x = x' = -L, (B7) 
27ri , else, 

as well as [<j)c(i~, x), <j> r .(T,x')\ = — [0c(tj %), 4>c{t, %')] — 27risgn(a; — x') , where sgn(0) = 0. Similar mode expansions 
were obtained in Refs. [36J. Given the mode expansion it is straightforward to calculate the four-point function (|F34|) . 
We find 

i/3i0c(2i) p iai^c(*i) p i«20c(z2) i/3 2 0c(z 2 )\ C(5 ai+a2:/3l+fe (zi - Z 2 ) 2aiQ2 (zi ~ Z2) 2/3lfe , , 

/UHP (z 1 -Z 1 ) 2Q ift(zi-Z 2 ) 2 ^/3i( Zl -f 2 )2a 1 /3 2 ( Z2 _ £2 )2a 2 /3 2 > I °J 

where C £ R is a constant which we set to one throughout this manuscript. This result implies Eq. (fT9|). The 
finite-temperature correlation functions are obtaine d 38 ! 58 by mapping (|B8|) onto a cylinder of circumference Vc/k^T. 



Appendix C: Calculation of the Green function: Spin sector 



The calculation of the correlation functions in the spin sector relies on the integrability of the sine-Gordon model 
on the half-line. We use the boundary state formalism introduced by Ghoshal and Zamolodchikov^ together with a 
form-factor expansion based on form factors obtained by Lukyanov and Zamolodchikov^. The analogous expansion 
for the quantum Ising model has been analyzed in Ref. (48j . We will first discuss the general formalism and then 
derive ([2%|) . 
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FIG. 13: Two-particle scattering and scattering off the boundary. 



1. Boundary state formalism and form- factor expansion 

Let us consider the sine-Gordon model ([14")) in the half-plane (r, x), r € R, — oo < x < 0. The boundary is located 
at x = and r denotes imaginary time (r = it). The Hilbert space of states associated with the semi- infinite line 
r = const, — oo < x < 0, is denoted by T-L^. We obtain the Euclidean action in its standard form by rescaling the 
hclds according to $ s — > $g = § S /K S and 6 S — > Gg = A' S S . The action of the sine-Gordon model with a boundary 
is then given by22 (we set v s = 1) 



SbsG = 



1 

loTT 



dr 



dx 



(d r <t>i s ) 2 + (d x <S>i) 2 -^cos(K s & s 



9b 



, K s 
I —I 



KL= f) -*° s /Ks 



(CI) 



where g s , g\, and $° arc free parameters. (We use the conventions < K s < 1, the action as given in Ref. [29[ is 
obtained by another rescaling of the hclds by \/8tt.) The cases <?b = and gh ~^ oo correspond to free and fixed 
boundary conditions, respectively. We stress that in the case of fixed boundary conditions $' s (x = 0) = &®/K s 
implies $ s (x = 0) = $g in the original system (TH|) . As was conjectured by Ghoshal and Zamolodchikm^ and shown 
independently^ by Maclntyre and Saleur, Skorik and Warner, the classical sine-Gordon model on the half-line (|C1[) 
possesses infinitely many integrals of motion and is hence integrable. 

We start by summarizing some results for the bulk sine-Gordon model, i.e. the theory without boundary. In the 
repulsive regime (K s > l/v2) a basis of the Hilbert space H is given by scattering states of solitons and antisolitons 



\0i 



zt 1 (6 1 )...zl n (e n )\o) 



• Al = 



■ • • Z ai (8i), 



(C2) 



where aj = ±1 and |0) is the ground state in absence of a boundary. Solitons and antisolitons are created by the 
operators Z_(ff) and Z\(9). They are characterized by a topological U(l) charge —1 and 1, respectively, while their 
energy and momentum are parametrized in terms of the rapidity 6 by E = A cosh 9 and P = A sinh 9. The dependence 
of the soliton mass A on the bare parameters in the action was obtained in Ref. 64[. We note that in the attractive 
regime (K s < 1/V2) breather (soliton-antisoliton) bound states occur as well. The operators Z a and Z\ fulfill the 
Faddecv-Zamolodchikov algebra^ (see Fig. [T5|) 



Za 1 {9i)Z a2 (92) 

Zl(9i)Zl(9 2 ) 
Z ai {9i)zU9 2 ) 



S b a \ b a l(9 1 -9 2 )Z b2 (9 2 )Z bl (9 1 ), 



- S b a \ b * 2 (9 1 ~9 2 )zl{9 2 )zl{9 1 ), 



(C3) 



2ir6{9 1 



^ aia2 +S b a ll]{9 



)zU9 2 )Z bl {9 1 ). 



Here S^\°^ {9) is the two-particle scattering matrix, which was derived in Refs. [32U66| . The unitarity condition reads 



S c a \Z(9)S b l 



bib 2 
<■■> ' 



S bl 5" 2 



Its non-vanishing elements are 



S++(9) = SZZ(9), S+Z(9) = Sim), Sil(9) = Sli(9), 



(C4) 



for which explicit expressions can be found for example in Ref. [44| . At the LEP (K s = 1/ V2) the scattering matrix 
simplifies to S blb a 2 2 (9) = —S b \S b2 2 , while in the spin symmetric case (K s = 1) one has^ 



S+U9) - S (9) 



r(i + g)r(§-g) 
r(i-f)r(i + f) 



S±Z(9) 



S (9), Sl7(9) 



9 - 



S (9). (C5) 



We note that the Faddeev-Zamolodchikov algebra (|C3[) is invariant under the unitary transformation Z a {9) — > 
e llp Z a (9), which changes the basis of scattering states. In terms of the basis states (|C2[) the resolution of the identity 
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reads 

id = 1°) <°I + E^E %^ • • • .fl>a„,..., 0l 0l -- a "(fl, • • ■ ,*n| . (C6) 

„=1 n ' {„,} j -» 

The boundary can be introduced^ as an infinitely heavy, impenetrable particle B sitting at x = 0. The ground 
state in presence of the boundary can then be represented as | Ob) = B |0). Scattering of elementary excitations off 
the boundary is encoded in the relations (see Fig. [13} 

zt{e)B = R b a (e)zl(-e)B, (C7) 

where the functions R b a {9) are the single-particle reflection amplitudes. In order to preserve intestability, the boundary 



reflection matrix R{9) has to satisfy a number of conditions which were discussed in Ref. |29[. At the LEP and for 



Dirichlet boundary conditions $ 8 (r, x — 0) = $g , < $!? < 7r in the original system ([T4]) . it is given by^ 

cosh( i~ ± i^f- 

= 7 sr-^ J^((?) = o. (C8) 

cosh (if ± i-f - 

For 7r/2 < $° the reflection amplitude i£f possesses a simple pole in the physical strip < 3m 9 < n/2, which 
indicates the existence of a boundary bound state. The overall sign of the reflection matrix is fixed by the requirement 
— iRes[i?±(0), 9 = ±i($° qp tt/2)] = -2cos$° > 0, see App. IC 51 below. Explicit representations of R for general K s 
can be found in Rcfs. 29.46.47]. For K s = 1 and $° = one finds in particular 

Rt(9) = - r ^f jT) ^ 2¥ ( cosh \ + isinh I) ■ = 0- ( C9 ) 

The vanishing of the off-diagonal amplitudes R±(9) = is a consequence of fixed boundary conditions and holds for 
general if s . 

Our aim is to calculate the time-ordered two-point function 

C{t, Xi ,x 2 ) - (0 b |r r Oi(r, a;i )O 2 (0,x 2 )|0 b ). (CIO) 

Here the time-dependence of the operators is given by Oi(r, x) = e rHh Oi(0, x) e~ rHb , where H\, is the Hamiltonian of 
the system in the presence of the boundary (|14[) . Given that in the Euclidean formalism r and x are interchangeable 
one may equally well designate x to be the "Euclidean time" . In this picture the equal-time section is the infinite line, 
x = const, — oo < r < oo, and the associated Hilbert space T-L is that of the corresponding bulk theory. The boundary 
at x = now appears as an initial condition which is expressed in terms of a "boundary state" |B). It was shown by 
Ghoshal and Zamolodchikov^ that the correlation function (|C10j) can be expressed as 

C(r,x u x 2 ) = e-* <°' % gill^ffi^ l B > . (C H) 

(0| a) 

Here s(Oi) denotes the Lorentz spin of the operator Oi, T x is the cc-ordering operator, which orders the largest Xi 
to the right, and |0) £ T-L is the ground state of the model on the infinite line. The spin-dependent phase factor is 
due to the rotation in Euclidean space; it was for example observed in the Green function of the Ising model with 
a boundary magnetic field^. As we have interchanged space and time and x is running from to — oo in the new 
framework, the r- and x-dependence of operators Oi(r, x) is now given by 

Oi(r, x) = e- xH e" irP O,(0, 0) e irP e xH , (C12) 

where H is the Hamiltonian of the system on the infinite line — oo < r < oo, and P is the total momentum. 
The boundary state, which encodes all informations on the boundary condition, is given by 

|B) = cxp Q jT ^-K a \0Zl(-0zi(0) |0) , (C13) 

where K ab {£) = i?|(i7r/2 — £). For example, the boundary reflection amplitudes K stated in ([2"Tj) and (|3"5|) are directly 
obtained from (|C8|) and (|C9[) . For general K s the amplitude K ah satisfies the boundary cross- unitarity condition^ 

K ab (Q = S a J(2H)K d %-0. (C14) 
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Furthermore, for fixed boundary conditions we have K^^(£) = 0. 

Below we calculate the spin part of the Green function (|Blj) using the boundary formalism presented above. 
Specifically we will evaluate the correlation function (jCllj) . where the operators are the soliton-creating and 
-annihilating operators e^^ 8 and e s 1 ^ 8 respectively. We define the n- particle form factor of an arbitrary operator 
O as 

fZ, ...,»„ :.,e n ) = {o\o \e u . . . A) 01 ,..., a „ = (o| o zt^) . . . z\ n {8 n ) |o) . (ci5) 

The form factors have to satisfy a set of relations, the so-called form-factor axiom s 39 i 41 i 43 , which we state for com- 
pleteness in App. IC 2l As the operators e~^^ s and e^ s create one soliton, their respective form factors (|C15[) vanish 
unless a,i — — 1. The form factors containing up to three particles were derived by Lukyanov and Zamolodchikov^ 2 -. 
In our conventions the single-particle form factors are given by 



(0|e-^ B \0)_ = sfz~ x e^ e 9 ' 4 , (0\e^\9)_ = ^e-^e- e/ \ (C16) 

where the normalization constant Z\ (not to be confused with the Faddeev-Zamolodchikov operators Z±{9) and 
Z±(6)) depends on K s and can be found in Ref. [Hf . Evaluation at the LEP yields Z\ w 3.32052 A 5 / 8 whereas at the 
SU(2) invariant point one finds Z\ sa 0.921862 A 1 / 2 . The three-particle form factors are known in terms of contour 
integrals, which can be explicitly evaluated at the LEP: 

(0|ei*-|9 1 ,9 2 ,9 3 )__ + J V2 cosh V ™* V* 

The three-particle form factors for other orderings of the U(l) indices can be easily obtained using the scattering 
axiom stated in App. IC 21 

The correlation functions to be calculated below contain matrix elements with incoming and outgoing particles, 

a ^-^{e x ,...Mo\u.---^i) bm ,... M , (cis) 

which possess kincmatical poles whenever 0i = £j and di = bj. These matrix elements can be decomposed into a 
"connected" and "disconnected" contributions. The latter are characterized by the appearance of terms like 5(6i 
signaling that some of the particles do not encounter the operator O in the process described by the matrix element. 
We deal with these terms following ideas by Smirnov 3 ^ that allow us to analytically continue form factors. Let 
A = {9i, . . . , 9 n } with 6i < 02 < ■ ■ ■ < n and *B = {£ m , . . . , £1} with £ m > £ TO _i > ... > £i denote two sets of 
ordered rapidities and introduce the notations 

Z[i] ai ,„ an = Z ai (6i)Z a2 (6 2 )...Z an (0 n ), (C19) 
&[% m ... bl = zl m {U)zt m _ 1 {U-i).--Zl i {ti)- (C20) 

Now let A\ and A 2 be a partition of A, i.e. A = A1UA2, where A\ contains n(A\) = n — k rapidities. As a consequence 
of the Faddeev-Zamolodchikov algebra we have 

z[i] ai ...a„ = s(l\Z) c a \;X z[Zu...o k z[Zt] Cfc+1 ... c „, (C21) 

where 5(^t|Ai) is the product of two-particle scattering matrices needed to rearrange the order of Faddeev- 
Zamolodchikov operators in Z[A] to arrive at Z\A 2 \Z\A^\. For example, if A = {#1, . . . , #4} and A\ = {6*2,6*3} 
it is given by 



Similarly we have 



Finally we define 



S(i\Zr-;;;X = 5 a\ (02 - fc) (63 - fc). (C22) 

Z\B] bm ... bl = Z^%] dm ... d!+1 Z^ 2 } dl ... dl S{%\B) d b Z:t- (C23) 

n 

S[A, B] ^. ;; »n = 5 nm J[ 2ir6 ajbj 6{8j - &■)■ ( C24 ) 
3=1 
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We are now in the position to analytically continue matrix elements as 

(o\z[t] ai ... an oz^] bm ... 6l |o) = sCA^y^x s{%\% d c:.t &\AX\% :: % 

A=Ai UA 2 
B = B 1 UB 2 

x(0\Z[i 1 +iQ] Ck+1 ... Cn OZ^^ 1 } dm ... dl+1 \0). (C25) 

Here the sum is over all possible ways to break the sets A and B into subsets and ~A*\ + iO means that all rapidities 
in A\ are slightly moved into the upper half-plane. Similarly, we could choose to analytically continue to the lower 
half-plane 

(o\z[i] ai ,„ an o&[% m ,„ bl \o) = s $&)l\:X 5[ZX] Z - Z 

B = B 1 UB 2 

x (0\Z[% - iO] Cfc+1 ... c „ O Z^r]^...^ |0). (C26) 

The factor d,A 2 (O) is due to a possible semi-locality of the operator O with respect to the fundamental fields creating 
the excitation o 39 i 41 ' 43 i 52 . If we use the operators defined in (|C29I) as fundamental fields and denote the mutual 
semi- locality factor of O and by l±(0), it is given by 

n(A) k k 

d A (0) = H l ai (0) => d A2 (e-^«) =n^, d A2 (eM°) =\[e^ a \ (C27) 

i— 1 i—l i—1 

The remaining matrix elements in (|C25[) and f|C26p can be evaluated using crossing 

(o|z[ti ± io] Cfc+1 ... c „ o zt[s 1 ] dm ... di+1 |o) = c ^— c "(^ +1 ±io, . . . , e in ±io\ o |e Jm , . . . ,& !+1 > dmj ... id(+1 

= dA 1 (0)C r c fc+1 e fc+I • ■ -C c „ e „ / e ° fc + 1 ,..., e „, dm ,..., d! + 1 (^ +1 +^±< + 1 ,- ■ - , ^„ + iTT ± i^„ , ^„ , . . (C28) 

where Cat = S a +b.o is the charge conjugation matrix and r/i — > 0+. The analytic continuation of general matrix 
elements (|C18|) with arbitrary orders of the rapidities can be obtained using the scattering axiom (see below) . 



2. Form-factor axioms 

For completeness we state here the used form-factor axioms. We follow Delfino^ 3 -. The n-particle form factor of an 
arbitrary operator O was defined in (|C15|) . We use the local bosonic fields 

0±(t,x) =exp(r^- s J dTd x & s (r,x)y (C29) 

as fundamental fields for the creation of solitons and antisolitons. The corresponding creation and annihilation 
operators are Z±(9) and Z±(6) introduced in (|C2[) . The form- factor axioms read: 

1. The form factors f~ a . . . ,9 n ) are meromorphic functions of 6 n in the physical strip < 3m9 n < 2tt. 
There exist only simple poles in this strip. 

2. Scattering axiom: 

/°,...,a i ,a i+ i,...,a„(01i • • • > • • • > On) 

= Sata£!i(0i ~ /oi,... 1 6 i +iA,-.o»(^l) ■ ■ ■ ■ ■ -,6 n ), 

with the scattering matrix Sala^iOi — Oi+i)- At the free-fermion point it is given by S^^(9) = —5 h a \5 a 2 2 . 

3. Periodicity axiom: 

/£,...,«„ (ft + 27ri, e 2 ,...,e n ) = i ai (o) fg (e 2 ,...,e n ,e 1 ), 

where l±(0) is the mutual semi-locality factor between the operator O and the fundamental fields . In 
particular, we have l±{e~^ B ) = =pi and l±(e^^ s ) = ±i. 
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4. Lorentz transformations: 

/£,...,«„ (01 + A, ■ • ■ ,0 n + A) = e s (°) A /£ B (0i, . . . , n ), 

where s(0) denotes the Lorentz spin of O. Here we have s(e ± 5*=) = 1/4 and s(e ± ^^ s ) = —1/4. 

5. Annihilation pole axiom: 

Res[/ a ° Ml) ... ia „ (0', Ml, ... , 9 n ), 6' = e + in] 

= i a, /£,...,„„ ... , 9 n ) [S b a \ . . . 5 b Z 5l - l a {0)S? b a \ (9 - e^Sttl (O-02)-.. S c Cn _± (9 - 9 n ) 

with the charge conjugation matrix C a t = <5a+fc,o- If there do not exist bound states in the model, i.e. for 
Ks — 1/2, these are the only poles of the form factors. 

We note that the precise form of the axioms depends on the basis of scattering states and thus changes under a unitary 
transformation of the operators Z a (9). 

3. Correlation functions 

In this appendix we derive (|20[) using the boundary state formalism. We start with the spin part of (|B1[) . After 
the rotation in Euclidean space this is given by (|C11|) . We insert a resolution of the identity (|C6|) and expand the 
boundary state (|C13[) in powers of K. This yields the double expansion (r > 0, x\ < x 2 ) 

00 00 

e -i/^s(0,x 2 )\ = (Q | g-i/^Cr.xO e -i/^ 8 (0,X3) | B) = 6aa , J- C n2m (T, Xl ,X2), (C30) 



I s 

n— m= 

where we have used s(e ± ^* s ) + s{e ± ^ B ) = 0. The operators e^^' and e^^ 3 change the U(l) charge by =Fl and ±1, 
respectively. As the boundary state has vanishing U(l) charge for Dirichlet boundary conditions (K aa (£) = 0) the 
correlation function is diagonal in spin space. Furthermore we have defined the auxiliary functions 

r u «. ^ 111/"°° <£i ■ ■ ■ d£ m f°° de x ... d6 n aibl(f , K a m b m(fz \ 

Cn2m{T ^ X2) -^^J^-^^L^^ K m- k {u) (C31) 



■i/^(r,xi) \ 0ni . . . j ft) «.™.c» (ft , . . . , 0„| e-kUU^,) 6; . . . ; _ u &„> 



,6„ 



We use the notations ^= +, ^= — , a = — for a = + and vice versa. We label the various terms in the expansion 
(|C30I) by the numbers of particles in the intermediate state n and in the boundary state 2m, respectively. The r- 
and x-dependence of the operators is given by (|C12|) . We have already assumed x\ < X2 to avoid additional phases 
due to the mutual semi-locality of the operators. For the calculation of the LDOS we have to take x\ — > x 2 — finally. 
The second matrix element possesses kinematical poles which we treat using (|C25[) . This introduces a third, finite 
summation in (|C30|) . which labels the "connectedness" of the corresponding terms. We note, however, that (|C'25|) 
and (|C26j) yield the same results. 

Let us start with the first non- vanishing term in the series (|C30|) , which is using (|C12[) given by (we recall that the 
center-of-mass coordinates are defined by R — {x\ + x 2 )/2 < and r — x\ — x 2 < 0) 

/OC J/J />00 lf\ 

— (0| e -V'+* \9) c c (9\ e -y°*« |0) e - rcoshe e iA ™» he = Z 1 ji / — e - rcoshe e iArsinhe . (C32) 
-00 27r J-00 27r 

We can rewrite this by shifting the contour of integration as 9 — » 9 + hr/2. The contributions of 9\z9 = ±00 vanish 
due to the exponential factors. As there are no poles in the strip < 3m 9 < ir/2 we find 

/OO if) y 
^ e i^,si„h9 e -Arcosh 9 = figif Ko (A^T 2 + r 2 /vl), (C33) 
-oo 27T 7T 

where Kq denotes the modified Bessel function^. 

The first term containing the boundary reflection amplitude K is C\ 2 . For / CT = — 1 it reads 

Cia = \ [°° ^r?~K ab {t) (0| e*^ \9) c c (9\e^ K,£) ab e^ rcoshe e 2 ^ 12 cosh « e iA " inh «, (C34) 
2 ./_„ 2tt 2tt 
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The first form factor vanishes for c ^ + and can be evaluated using (|C28|) 

(0| e* 03 \0) + =+ (0\ e'^' |0)* = e 1 ? (0| e"^= |6> + Itt)* = \fz~ x e^ e e/i . (C35) 
For the second matrix element we use (|C25[) . which explicitly yields 

+ (6\ e± A K,0-+ = + (* + i0| e**" |-f,0-+ + 2 ^(0 - (0| e^ s K>- + 2 ^ + £±^(-20 <0| e* A , 
+ <0| e**- = + + i0| e* A + 2^(0 + S+Z(-20 (0| e**- |£>_ ■ 

(C36) 

This leads to two contributions which we denote by C® 2 and C| 2 respectively. The additional upper index denotes the 
number of lines connecting the operators (the "connectedness" ) , i.e. the number of internal ^-integrations left after 
using (|C25p . The first terms on the right-hand side of (|C36J) in each equation together yield C\ 2 . We will calculate 
this term at the LEP in the next section. On the other hand, the disconnected piece is given by 



,2^ 



K-+{0 + K+-{-i) Stl(2£) + K~ + (-0 St+(2£) e«/ 2 e 2 ^ Rcoshi e iA ™» h «. (C37) 



Using the boundary cross- unitarity (|C14[) the terms in the square brackets equal 2K + (£). With the similar calculation 
for f a = l we arrive at 

C? a = Z X r ^ K™(0 e«/ 2 e 2 ^ Rcosh t e iATSi ^, (C38) 
J-oo 2tt 

The second term in (|20p is now obtained by shifting the contour of integration as £ — > £ + i7r/2 while noting that for 
the boundary condition $ s (x = 0) = the reflection amplitude does not depend on a and is analytic in the physical 
strip < 3m£ < tt/2. If $> s (x = 0) 7^ 0, however, the reflection amplitude may possess a pole in the physical strip. 
We will calculate the resulting term in App. IC 51 



4. Higher-order terms 

In order to estimate the truncation error in (|20[) , we will calculate the leading corrections due to a higher number of 
particles in the intermediate state as well as higher-order corrections due to the boundary. The resulting corrections 
to the LDOS arc discussed in Sec. IIV CI We will restrict ourselves to the LEP, where the form factors are given by 
CHI, S%%{9) = and K<*(t) = -K ba (-t). 

The leading correction due to a higher number of particles in the intermediate state is given by C30, 

° 3 ° ~ Zl ~ "T^F" cosh 2 E=S C osh 2 S=S 6 " s 6 ^ ■ (C39) 

The resulting contribution to the LDOS discussed in Sec. IIV Cl is denoted by N30. We note that C20 = 0. 

The first sub-leading term due to the boundary is given by C\ 2 , i.e. the connected piece of C12 obtained from the 
first terms in (|C36|) . For f a = —l this term reads using (|C35|) and (|C17[) 



r° djdd K + (£) smh e+ Z+™ A rcosh e 2A l2COS h{ iArsinhe rnAri\ 

Cl2 Zl J-00 ^ ^Shf COSh^±^ 6 ^ 6 ' (C40) 



where 77 — ^ 04-. We can handle the singularity at 9 = £ — irj by shifting — > + i?r/2. Performing the same steps for 
/cr = 1 we arrive at 



G\ 2 =Z x e -^ r ^^^ll_zi e ^ e ^^e^^ e - ATCOshe . (C41) 
12 V2 7_oo27r27r cosh£ ie« + e e v y 

The next term in the series (|C30[) is C32, its disconnected piece is similar to C\ 2 and reads explicitly 
,i . _ 7 e~' f f°° d£_ d8 K a,y (C) ie^ + e e s/2 iA rsin he „2A Xl COB he o -Arc OS he 



32 V2 i_ 0O 2 7 r27r cosh£ e«+ e - i v ; 



22 



The terms C\ 2 and C\ 2 are of the same order and yield together the contribution to the LDOS denoted by N\\. 
The term resulting in N 2 \ is 

n 2 _ if f°° d£ d9 1 d9 2 jA r ^. sinhflj „2iAj?, sin hg ^-Ar(T. . cosh S.+cosh Q 

"lv CT5 (C + if)e«/ 2 sinh^^i s i n h^- i A^ CT (£ + if ) er^ 2 sinh 2 ^ 



cosh 2 cosh ^ cosh ^ 2 [] . CO sh ^ cosh 



(C43) 



We note that after analytic continuation r — > ft + S and Fourier transformation t E the exponential factor 
g-Ar^cosh^+coshC) resu its in a vanishing of N 21 for energies A < 3A. 

The final term we wish to evaluate explicitly is the disconnected piece of C14. Considering first f a = — 1 and keeping 
in mind that we have restricted ourselves to the LEP, we can start with 

2 J-oo (27rj 2 2tt 

x(O|e**-|0> + + (0|e^=|-6,a,-6,6>- + -+- (C44) 
In the second matrix element we keep only the disconnected piece 

+ <0|ei^|-a,&, -&,&>_+-+ = 2n8(0-ti) (0| e** -£ a ,&>__ + + 27r<J(0-&) (0| e**- 6, -&>_+_ + • • • 

(C45) 

In the resulting term C° 4 we can shift the contour of integration, £1 — > £1 + i7r/2, to obtain 

r° --^L P -ii r K as (t. I \*\ K™(M eiil+i2)/2 e^+ie^ 2i A fisinh6 2 A X2CQsh g 2 A r corf.fi c C4fil 

° 14 " 72 * 1+1 *) K (6) cQsh6 e e . (C46) 

In the last step we have assumed that there exist no boundary bound states (see below). Finally, we mention that 
the next term, C34, equals C44 with the coordinates x\ and 2; 2 interchanged. These two terms together yield Nq 2 . 

The remaining two terms N12 and N03 discussed in Sec. IIVCI follow from C\± + C34 + C54 and C° 6 + C° 6 + C$ 6 
respectively 



5. Boundary bound states 

As discussed above, general Dirichlet boundary conditions $ s (0) = $g 7^ can result in the appearance of a 
boundary bound state. If A 2 7r < $^ < it the boundary reflection amplitude A h (£) has a pole in the physical strip 
< 3m£ < tt/2 located ai4^ £ = i(?r - $°)/(2 - 2A' 2 ). On the other hand, K+~(\) is analytic in the physical strip 
but has a pole for — n/2 < 3m £ < 0. We write the respective residues as 



iRcs 



K* ± (ZU = ±i<y]=B>0, 7=^-^, (C47) 



where B depends on A s only. We have checked the sign of B by performing an explicit mode expansion at the LEP 
as well as studying the spectral function of the correlator (0b| e~ la *^ T '^ e 10 *^ '^ |0b) (for which the relevant form 
factors were obtained in Ref. (o7j)- 

In the presence of a boundary bound state the poles of A T± (£) will contribute whenever we shift the contour of 
integration £ — > £±iV/2 in a given term in the form- factor expansion (|C30|) . The leading term of this type is obtained 
from (|C38|) , which yields (|41[) by a straightforward calculation. The sub-leading term can be obtained similarly from 
Cj 4 and C§ 4 . At the LEP it is given by (B = -2cos$°, tt/2 < $s < n ) 



U s -^)5^Z 1 VSe-^cos^ s e- 2 ^ Rcos ^e- AT ^ [°° ^ — >H - V e^ 2 e 2 ^ cosh ^. (C48) 
v y 2tt cosh£ sinh / £ _ i^ol 



cosh 


3 + 




sinh 
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FIG. 14: Branch cut and deformation of the contour of integration used for k < and k > respectively. 



Appendix D: Fourier transformation of the LDOS 



We calculate the auxiliary function 

r-0 



I(u>,k) 



dR 



dR 



dt- 



dt 



,\{uit-kR) 



1 



(v c t - 2iR) („ cT + 2LR) 



2R 

VrT 



2c 



T->it+(5 



(v c t-2R-i5) a ( Vct + 2R-iS) b W 



2R 



-id 



2c 



(Dl) 
(D2) 



where v c , a, b, c € M, v c > 0, a + b < 2 and c > —1/2. We substitute R — > —R and t — > — t, introduce s = v c t/2R and 
?7 — > 0+, and perform the resulting i?-integral (3.381.4 in Ref. (63|), which yields 



I(w,fc) 



^iTrfa+b— c) 



r(2-a-6) 



f/.S 



'^s - fc - in" 



a+b-2 



(s - 1 + i<S) a (s + 1 + W) 6 (s + W) 



2c ' 



(D3) 



For cj < the integrand has all its branch points in the lower half plane and the integral over s vanishes as long as 
c > —1/2. Hence we find (so = v c k/2co) 



I(u,k) = -Q(u)T(2-a-b)- 



(a+b-c) ^a+b-2 poo 



2 V, 



a+b-1 



ds 



(s - s - 177) 



a+b-2 



(s - 1 + iS) a (s + 1 + iS) b (s + iS) 2c ' 



(D4) 



First consider the case k < 0. The numerator of the integrand has a branch point at s = sq + irj in the upper half 
plane. We place the cut running from —00 + it] to so + if] with constant imaginary part (see Fig. I14[) . Now we deform 
the contour of integration and rewrite the integration above the cut as integration below the cut using 



ds (s — Sq — irj) 



a+b-2 



9(s) 



ds (s — so — irj) 



a+b-2 



1 - e 



2m(a+b) 



Assuming 1 < a + b and substituting s = so/i this yields for the integral in (|D4 



2i sin(7r(a + b)) \sq\ 



0+6-1 



dt 



t 2c (1 - 1) 



a+b-2 



(s - t + iS) a (s + t + iS) b (s + iS) 



2c ' 



Filially, using (recall so < 0, 5 = 0+) 

(s - t + iSy a = (s (l - (l/s + i5)t)y a = |s |- a e-™ (l - (l/s + W)tp 
(s +t + iSy b = {s (l + (l/s -i5)t)y b = \s \- b e- M {l + (l/s -iS)ty b 



-2c 



(so+iS) = 1 s 1 



2c „-2iric 



as well as e 



27ric /|s | 2c+1 



-(l/so — iS) 2c+1 and r(z)r(l — z) = it/ sin(7rz) we obtain 



I(u,k< 0) = 



^e(w)e- i t( 2c ~ 1 ) io a+b - 2 ( 1 



r(a + 6-l) u^ 6-1 \s 



iS 



2c+l .1 



t 2c (1 - t) 



a+b-2 



(l-(l/s + i<S)i) a (l + (l/s -i<5)t)' 



For fc > we place the cut as shown in Fig. [14] Performing the same steps as above we find 



I(u,k> 0) = 



r(a + 6-l) w^- 1 V s o 



— +w 



2c+l .1 



dt 



t 2c (1 - t) 



a+b-2 



(l-(l/s -i<S)i) a (l + (l/s + i<S)t)' 



(D5) 



(D6) 



(D7) 
(D8) 
(D9) 



(D10) 



(Dll) 
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Wc can write (ID10[) and (|D11|) together as 

I(w, fc) = * Q ^ )e r( ^T h lo^ C + ^ " 2C+1 Ji (2c + 1, a, M + 6 + 2c; u*, -u) , u = ^ + i sgn(fc) 8. (D12) 

1 (a + o + 2c) «c 

Here we have used the integral representation (|E2[) of Appell's hypergeometric functions 9 -, which is valid for 1 < a + b. 
Analytic continuation in the parameters a, b, and c then yields I(u), k) for a + b < 2 and c > —1/2. At K c = 1 one 
finds F 1 (2c+l,a,b,a + b + 2c;u*,-u) = F 1 (l, 1/2,0, 1/2; it*, -it) = 1/(1 - «*). 
In the same way one can show 

/-oo gi(wt-fc-R) (_j^a+&+2c (2R) 2c 

dR I dt 



(v c t -2R- i5) c (v c t + 2R- iSf (v c t - i5) a + b 



as well as (A > 0) 



[°° e iM-fcfl) (_ i) a+6+2c e Afl / 2i? X2c ^Q^) g-if (2c-l) F(2 C + 1) 



t-2R-iS) a (y c t + 2R-iS) b W-iSJ ^—"v^- 1 T(a + b + 2c) 

x (^+isgn(*)*) ^ F D ) ( 2c + 1 '«^' 2c - a + & + 2c ;ai- i T I -a- i T ; - i T) ; 

where denotes Lauricella's hypergeometric function of three variables (see App. |E]) . For a = 1/2 and 6 = c = 
(|D14[1 simplifies to 2\tJwv^Q{lo) / yjuj / (v c k — 2lu + iv c A). 

Appendix E: Hypergeometric function of several variables 

Hypergeometric series of several variables were first studied by Lauricella^. They are defined by 

F£\a,p u ...,Pn,r,z 1 ,...,z n )= jr (Q)mi+ -;"!" (/3l)mi --- (/3 " )m " z ™\--- z ™ n \ Zi \<l. (El) 

mi ,.tL=o (7)m 1+ ...+m„ mi!---m„! 

The special cases^ 9 - n = 1 and n = 2 are Gauss hypergeometric function F^ = F(a, ft 7 ;z), and Appell's hyper- 
geometric function F^ = F±(a, ft, ft, 7 j z i, z i), respectively. The function F^ possesses the Euler-type integral 
representatio n 49 ! 68 

A"), a o .... . N_ r( 7 ) f\ u p-l(l-^)T-«-l 



(a, ft,..., /3„, 7 ;z 1 ,...,z„) = — — ^ -/ dt- ^ -rg-, <Kea>0, SRe( 7 -a) > 0. (E2) 

r(a) r( 7 - a) Jo (1 - zii) Pl ■ ■ ■ (1 - z„t)P™ 

Furthermore the following relations hol d 49 ' 56 

f£° (a, ft,..., ^ 7 ; Zll ... A ) = (l-^.-(l-^< ) ( 7 -a, ft,..., ft, 7; . . . , (E3) 

Fi(a,/3i,/3 2 , 7 ;l,l)= 7; n g ^ for 7 ^ 0, -1, -2, . . . and 7 > a + ft + ft. (E4) 

F( 7 -a)r( 7 -ft - ft) 

Appendix F: Properties of N>(E, 2k F + q) 



In order to analyze the dispersing features and singularities of (|26|) we first note that F\ (2c+l, a, b, a+b+2c; u*, —it J 
possesses singularities at it = ±1. 

Let us first study Nf . The integrand has singularities at 

(i) E- A cosh 6 = 0, (ii) 2(£ - A cosh 0) = ±v c q. (Fl) 
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Inserting (i) into (ii) immediately yields a feature at q = 0. Using (|E3|) and (|E4[) one can extract the g-dependence 
{v c q) a+b ' 2c l to obtain (f29|) . On the other hand, (i) will be stationary at 9 w 0. Inserting this into (ii) directly yields 
the dispersion relation ([30j) . The suppression of the dispersing peak for g < follows from the relative strength of the 
singularities at u = ±1. 

In the same way the integrand in A^" has singularities at 

(i) E - A cosh 6» = 0, (hi) 2v s (E - Acoshtf) = ±v c (v s q - 2Asinh0). (F2) 

Inserting (i) into (hi) directly yields the dispersion relation (|31[) . Furthermore, we can rewrite (hi) as 

(iv) cosh T ^ S inh 9. (F3) 

A 2A t> s v ; 

If and only if v c < v a , the right-hand side in (iv) becomes stationary at 9 = 9 = ±areosh(u s / yuf — wf) ■ In principle, 
this leads to the relation (|32|) for arbitrary g. However, this dispersing feature only exists when — arcosh(^) < 9 < 
arcosh(^). Together with (j3"2"j) this yields the condition q < \q\. 

Finally, in order to prove (|45|) we use that (2c + 1, a, b, 2c, a + b + 2c; w|, — M3, — u'^j is regular as E — > i?bbs+, 
which directly yields a = l — a — b — 2c=l — l/(2K%). 
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